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ABSTRACT 


Let  denote  the  Lucas  numbers,  i.e.  the  sequence  of  integers 
L^  satisfying  the  recurrence  relation  L^^^  = L^  + L^  p with  L^  = 2, 

I*!  = I-  From  various  sources  the  following  conjecture  was  formulated; 


The  number  n is  a prime  if  and  only  if  (1)  L^  = 1 (mod  n). 

In  the  reference  [3]  Hoggatt  and  Bicknell  have  shown  that  the  "only 
if”  part  is  correct,  while  the  "if”  is  wrong,  counter  examples  being  the 
numbers  705,  2465,  2737  and  others,  all  of  which  are  composite  and 
satisfy  (1).  In  this  paper  we  derive  these  results  of  [3]  making  extensive 
use  of  the  matrix  approach  to  Fibonacci  numbers  as  described  in  the  book 
[2,  Chap.  11].  We  also  describe  the  results  of  extensive  computations  due 


to  George  Logothetis  and  done  in  November  1976, 

AMS(MOS)  Subject  Classification  - 10A35 
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Work  Unit  No:  2 - Other  Mathematical  Methods. 
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ON  THE  MATRIX  APPROACH  TO  FIBONACCI  NUMBERS 
AND  THE  FIBONACCI  PSEUDOPRIMES 

Jack  M.  Pollin  and  I.  J.  Schoenberg 

Introductiot^  We  consider  sequences  (x^)  of  Integers  satisfying  for  all  n the  re- 
currence relation 

(1)  = X + X , . 

n+1  n n-1 

The  x^  are  uniquely  defined  if  we  prescribe  the  elements  of  the  "initial  vector"  (x^jX^).  On 

choosing  (Xq,Xj)  = (0,1)  we  obtain  the  Fibonacci  numbers  while  the  choice  (x^,Xj)  = 

(2,1)  gives  the  Lucas  numbers  x = L . 

n n 

In  [3]  V.  E.  Hoggatt  and  Marjorie  Bicknell  discuss  the  following  conjecture  of  K.  W. 
Leonard  (unpublished). 

Conjecture  1.  We  have  the  congruence 


(2)  L^Hl(modn),  (n  > 1) 
if  and  only  if  n is  a prime  number. 

Among  the  many  interesting  results  of  [3]  we  single  out  the  following. 

Theorem  1.  The  "if"  part  of  Conjecture  i is  correct,  l.e. 

(3)  ' ^p  ^ ^ (mod  p),  where  p is  a prime. 

Theorem  2.  The  "only  if"  part  of  Conjecture  1 is  wrong,  as  shown  by  the  congruence 


(4)  = 1 (mod  705)  , 

while  705  = 3.  5.  47  is  composite. 

We  owe  to  D.  H.  Lehmer  an  informative  letter  [4]  in  which  he  expresses  familiarity  with 
these  results;  also  that  composite  numbers  n that  satisfy  (2)  are  called  Fibonacci  pseudoprimes, 


which  we  abbreviate  to  F.  Psps.  In  [3]  the  authors  report  on  the  basis  of  computer  results  that 
beyond  705  the  next  F.  Psps  are 

(5)  2465,  2737,  3745,  4181. 

Conjecture  1 was  communicated  to  one  of  us  several  years  ago  by  Richard  S.  Field,  of 


Los  Angeles.  We  became  aware  of  the  paper  [3]  only  recently.  Before  this,  in  November  197  6, 
Sponsored  by  the  United  States  Army  under  Contract  No.  DAAG29-7  5-C-0024. 


Geors®  Logothetis,  a graduate  student  in  Computer  Science  In  Madison,  using  Professor  George 
Collins'  SAC  2 program,  found  for  us  not  only  the  five  F.  Psps  already  mentioned,  but  also  the 
two  new  ones 

(6)  5777,  6721, 

and  that  these  seven  numbers  are  the  only  F.  Psps  which  are  ^9161, 

In  the  present  paper  we  do  the  following. 

1.  We  present  a proof  of  Theorem  1 that  uses  from  elementary  number  theory  only  Euclid's 

Lemma. 

2.  We  give  a second  proof  of  Theorem  2 and  also  establish 

Theorem  3.  ^z4bS  ^ ^ ‘ 

These  numerical  results  are  here  derived  by  the  matrix  approach  as  described  in  (2, 

Chapter  11].  In  [3,  p.  211]  Theorem  2 is  proved  in  a few  lines  by  showing  that  the  sequence 

L mod  705  has  the  period  704.  Since  L,  = 1,  the  relation  (4)  follows.  In  §3  we  describe  this 
n i 

method  of  periods  and  show  that  while  it  proved  Theorem  2,  it  does  not  work  to  establish  Theorem 


3.  In  [4]  D.  H.  Lehmer  stated  that 

(7)  2737  = 7.  17.  23  is  a Fibonacci  pseudoprime, 

and  that  the  method  of  periods  will  apply.  This  we  verify. 

3.  In  §5  we  show  that  the  matrix  approach  allows  us  to  develop  ab  Initio  some  of  the 
basic  properties  of  Fibonacci  numbers  as  presented  in  [1,  §10.14],  As  we  assume  no  previous 
knowledge  of  Fibonacci  numbers,  this  paper  may  serve  as  an  introduction  to  these  numbers. 


4,  The  failure  of  the  "only  if"  part  of  Conjecture  1 suggests  a search  for  classes  of  com- 
posite numbers  n which  are  not  Fibonacci  pseudoprimes.  In  §6  we  state  some  modest  results 
in  this  direction.  These  suggested  the  following 
Conjectur^^.  _If_  n > 1 , then 

(8)  L 1 (mod  n^)  . 

n 

Again  G.  LogotheUs  showed  (8)  to  hold  for  n <7611.  Some  further  striking  results 
obtained  in  the  course  of  this  computation  are  described  at  the  end  of  the  paper. 
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1.  A proof  of  Theorem  1.  Observe  that  the  Lucas  numbers  L are  explicitly  given  by 

(1.1)  = {--fi- ) + (■^^) 

because  (1  ± VT)/2  are  the  roots  of  the  characteristic  equation  - x • 1 = 0 of  (1),  hence 
the  right  side  of  (1.1)  satisfies  (1),  while  it  assumes  the  same  initial  values  as  for  n = 0 
and  «n  = 1.  Let  now  n = p be  a prime  > 2.  Expanding  the  binomials  and  canceling  the  irra- 


tional terms  we  find  that 


2**  - 2 


Applying  in  the  numerator  of  the  last  term  the  binomial  expansion  of  (1  + 1)  , we  obtain 

The  left  side  is  an  integer,  while  the  right  side  is  of  the  form  p a/b  , where  p does  not  divide 
b,  and  therefore  (p,b)  = 1 . By  Euclid's  Lemma  we  conclude  that  b divide  a,  which  proves  (3). 
2.  The  matrix  approach  and  a proof  of  Theorem  2.  We  replace  the  relation  (1)  by  the 


vector  recurrence  relation 


X 0 1 X . 

(x"  ) ■(, . ;■*) 

n+l  n 


to  which  it  is  visibly  equivalent.  Writing 


= (”  M, 

' 1 1 ' 


and  iterating  (2.1)  we  obtain  that 


’‘ax 


This  brings  to  bear  on  our  problem  the  powerful  tool  of  matrix  multiplication.  To  prove  Theorem  2 


it  suffices  to  work  modulo  7 05.  We  observe  that  (2.  3)  implies 


(l  ^ f) 

' **705  ' ^ 


and  that  we  are  to  determine  the  matrix  A (mod  705).  This  Is  readily  done  with  a hand  calcu- 
lator If  we  use  the  binary  representation  of  704  which  Is 
(2.  5)  704  = 64  + 128  + 512  = 2^  + 2^  + 2^  . 

By  successively  squaring  of  matrices,  and  working  mod  705  throughout,  we  find  the  matrices 
2''* 

A <n»od  705)  for  Ic  = 1, 2, . . , ,9  , and  In-partlcular 

.2^  /142  423  \ .2^  _ / 283  141  \ .2’_/424  564  \ 

^ "(423565)’  ^ =(141424)'^  =(564  283}’ 

Multiplying  together  these  three  matrices,  mod  705,  we  find  by  (2.  5)  that 


J04  ^ / 142  423 

^ " ( 423  565 


^ , (mod  705)  . 


Now  (2.  4)  shows  that 


/ So4\  _/142  423v2\  /707  \,/2\ 

V "(423  565  A 1 ) (l411  ) Vl)»<"'°'^  • 

Therefore  - 1 (mod  705)  and  Theorem  2 Is  established. 

A few  remarks  on  these  matrix  operations  are  In  order.  Observe  that  A is  a symmetric 
T 

matrix,  l.e.  A = A.  We  also  know  that  the  product  BC  of  two  symmetric  matrices  that  com- 
mute (BC  = CB),  is  also  symmetric.  Since  any  two  powers  a"*  and  a"  clearly  commute,  it 
follows  that  all  powers  A*"  are  symmetric.  This  means  that  in  multiplying  two  powers  of  A 
we  need  to  compute  only  one  of  the  two  elements  off  the  main  diagonal. 

The  matrix  multiplications  performed  above  require  the  following  important  check  against 
errors.  Passing  to  determinants,  from  |A|  = -1,  we  conclude  that  Ia"*]  =(-!)*".  Since  above 
all  our  exponents  m are  even,  we  see  that  1a"*|  = 1,  and,  of  course,  also  1a"'|  = 1 (mod  705). 
The  check  is  to  verify  that  after  each  matrix  multiplication  the  resulting  product  M satisfies 


|M|  s 1 (mod  705). 


I 


3.  On  the  Hoggatt-Bicknell  proof  of  Theorem  2.  In  order  to  make  this  paper  self-suf- 
ficient we  establish  the  known  Lemmas  below.  Let  k be  given,  k >1,  and  let  us  denote  by 

(L  , mod  k)  the  sequence  (L  ) of  Lucas  numbers  reduced  mod  k . 
n n 

Lemma  1.  The  sequence  (L^>  mod  k)  is  periodic. 

Proof;  Clearly  (L^,  mod  k)  is  periodic  If  and  only  if  for  some  r and  s we  have 

**r’  ’‘r+l^  “ ^^s’  *s+P  r < s . 

It  follows  that  there  is  no  periodicity  if  and  only  if 

for  every  pair  (r,  s),  such  that  r < s we  have  (x  ,x  ) (x  ,x  ) (mod  k). 

r r 1 1 s s I X 

But  this  is  obviously  impossible,  as  there  are  only  k^  distinct  pairs  (u,v)  (mod  k)  available. 

The  Hoggatt-Bicknell  proof  of  Theorem  2 is  based  on  the  following  sufficient  conditions 

for  (L^,  mod  k)  to  have  the  period  m . 

Lemma  2,  If  the  following  conditions  are  satisfied 
t 


^ ""A  1 li  * j 

A is  a period  of  (L  , mod  a ) , 


Ajm  for  all  1 , 


(3.4)  m is  a period  of  (L^,  mod  k). 

Proof:  By  (3.  2)  L . . s L (mod  a.)  for  all  n . Bv  (3.  3)  it  follows  that 

n+Aj  n 1 

(3.5)  L , = L (mod  a.)  for  n , and  all  i , 

n+m  n 1 

because  a multiple  of  a period  is  also  a period.  Now  (3.1)  and  (3.  5)  imply  that  L = L (mod  k) 

n+m  n 

for  all  n , which  proves  (3.4). 

Lemma  2 applied  nicely  to  the  case  of  k = 705  = 3.  5.47,  for  (3.1)  holds  with  t = 3, 

3)  = 3,  a^  = 5,  a^  = 47,  Simple  direct  calculations  with  L^  show  that  (3.  2)  is  satisfied  with 

Aj  = 8,  A^  = 4,  Aj  = 32.  Also  (3.  3)  holds  for  m = 704  because  8,  4,  and  32,  are  all  divisors 

of  704.  By  Lemma  2 we  conclude  that  L s l (mod  705)  for  all  n . In  particular  for  n=l 

n+704  n 

we  obtain  = 1 (mod  705),  which  proves  Theorem  2,  For  n = 0 we  also  obtain  that  = 

Lq  = 2 (mod  705).  which  we  already  know  from  (2.7). 


I 


This  method  will  not  allow  us  to  prove  Theorem  3,  Indeed,  the  relation  (4.  3)  below  shows 
that  m = 2464  Is  not  a period  of  (L  , mod  2465). 


(4.1) 


4.  _Ajgr^f_o|_Th^rem^  By  (2.  3)  we  are  to  determine 

(mod  2465)  . 


From  2464  = 32  + 128  + 256  + 2048  = 2^  + 2^  + 2®  + 2^*  we  obtain 

ij.  *2464  .2^  2^  2®  2^^ 

(4.2)  A =A-A-A'A 

By  successive  squaring  of  matrices  mod  2465  we  find  that 


, 379 

1714 

a 

1393 

1886 

= Vni4 

2093/ ' 

A = 

(1886 

814  ) ’ 

( 495 

1482  \ 

2“ 

/ 1858 

1221  \ 

fM 

00 

III 

1977  • 

pr  = 

^ 1221 

614  • 

(4.3) 


( L.,,,  ^ ( 783  900  X 1 ) ■ (z466)  ( 1 ) (mod  2465). 


(mod  2465) 

A bC»  « «./  A ~X 

Multiplying  these  together  we  find  by  (4.2)  that 

117  783 

2464  _ / N 

^ “ W83  900  ) 

and  finally,  by  (2.  3) 

L 117  783  .,2  1017.  1017 

i no-i  ortA  K t J “ f 

2465 

Thus  L2465  ^ ^ (mod  2465),  which  proves  Theorem  3. 

The  information  that  = 1017  (mod  2465)  shows  that  m = 2464  Is  not  a period  of 

(L|^  , mod  2465),  and  this  is  the  reason  why  the  method  of  §3  would  not  work. 

Similarly  we  can  work  out  on  a hand-calculator,  such  as  SR-51A,  the  matrix  A^'Vmod  n) 
for  any  n < 10^.  Indeed,  all  matrix  multiplications,  mod  n,  are  leasable  because  all  numbers 
that  we  encounter  are  < 10^®,  the  capacity  of  the  calculator. 

In  [4]  D.  H.  Lehmer  pointed  out  that  the  second  number  of  (5),  hence  27  37  = 7. 17.  23 
is  a Fibonacci  pseudoprime  and  that  Lemma  2 applies  to  show  it.  This  we  easily  verify:  Lemma 
2 applies  to  k = 27  37,  with  t = 3,  a^  = 7,  a^  = 17,  a^  = 23,  Aj  = 16,  Aj,  = 36,  A^  = 48,  and 
m = 27  36.  Therefore  27  36  is  a period  of  (L  , mod  27  37)  and  it  follows  that  L , = 2, 

^27  37  ^ * (mod  27  37).  Therefore  (7)  is  established. 


5.  Further  applications  of  the  matrix  approach.  Our  applications  in  §§2  and  4 were 
mainly  computational.  We  now  wish  to  show  how  the  matrix  A allows  us  to  develop  ab  initio 
some  of  the  best  known  properties  of  the  Fibonacci  numbers. 

Let  us  make  the  relation  (2.  3)  or 


(5.1) 

more  explicit  by  writing 


II 

> 

o’ 

^n+1^ 

a b 


(5.2)  J>) 

n n 

whereby  it  becomes 


(5.3) 

This  easily  generalizes  to 

(5.4) 

Indeed  by  (5.1) 


X = a X-  + b X, 
n n 0 n 1 

X . , = c x„  + d X, 
n+1  n 0 n 1 


X = a X,  + b X,  , 
n+k  n k n k+1 

^ . 1. . 1 = c X,  + d X,  , 

n+k+1  n k n k+1 


i- 


( ) = a"+'^(  = a"  . a''(  = a"(  ) , 


n+k+1 


k+1 


again  by  (5.1).  Now  this  and  (5.  2)  show  that  (5.4)  holds.  We  obtain  x = F if  we  choose 

n n 

Xp  = Fjj  = 0 and  x^  = F^  = 1 and  (5.3)  shows  that 


(5.5) 


F = b 
n n 

F = d . 
n+1  n 


Applying  (5.4)  to  x^  = F^  and  k = 1,  observing  that  Fj  = 1,  = 1,  we  obtain 


F ,,  = a + b 
n+1  n n 

F . _ = c + d 
n+2  n n 


These  relations  and  (5.  5)  show  that 


®n  “ ^n+1  " ^n  ~ ^n-1  ’ 
°n  " “^n+2  ■ ^n+1  " ^n  ’ 


-7- 


We  have  thus  shown  that 


F F 

a"  = ( "-*  " ) . 


See  also  [2,  Theorem  II]. 


Our  previous  remark  that  |A  | = (-1)  shows  that 

which  is  a known  relation  derived  in  the  same  way  in  [2,  Theorem  III].  From  (5.6)  we  also  see 

that  the  elements  of  all  the  matrices  of  §§2  and  4 are  appropriate  Fibonacci  numbers  reduced  by 

the  moduli  705  and  2465,  respectively. 

Let  us  derive  the  known  property  that 

(5.8)  F divides  F if  r>0. 

n nr  — 

From  (5.4)  and  (5.6)  we  obtain  for  x = F the  relation 

n n 

(5.91  ( *^"  )(  "^^  )• 

'’k*! 

Replacing  here  n and  k by  nr  and  n , respectively,  we  obtain 


F F F 

n(r+l)  \ ( nr-1  nr 


X ; ) . 


whence 


^n(r+l)+l  ^nr  ^nr+1  ^n+1 


I"  / .,5  = F , F + F F . 
n(r+l)  nr-1  n nr  n+1 


This  shows  that  if  F divides  F , then  F also  divides  F , .,,,  and  this  proves  (5,8)  by 
n nr'  n n(r+l)’  \ i i 

induction,  since  (5.8)  is  obvious  if  r = 1 . 

As  a further  example  let  us  establish  the  known  property; 

(5.10)  If  (m,n)  = d then  (F  ,F  ) = F . . 

— m n u 

Since  d divides  m and  also  n , it  follows  from  (5.8)  that 


(5.11)  F.  divides  F and  also  F . 

dm  n 

There  remains  to  show  that  F,  is  the  greatest  c.d.  of  F and  F . Let  r and  s be  such 

d •* m n 

that  d = mr  + ns.  From  (5.9),  on  replacing  n and  k by  mr  and  ns,  respectively,  we  obtain 


-8- 


i'. 


mr+ns 


)■  ( 


f 1 I’  s 
mr-1  mr  N 


mr+ns+1 


F F t 

mr  mr+1  ‘ns+1 


) { r ) 


This  shows  in  particular  that  F , = F can  be  written  as 

d mr+ns 


(5.12) 


F = F F + F F 

d mr-1  ns  mr  ns+1 


By  (5.8),  any  divisor  6 of  F and  of  F also  divides  F and  F , and  by  (5. 12)  that  6 

m n mr  ns 


also  divides  F,.  Therefore  F,  is  the  greatest  c.  d.  of  F , F , and  (5.10)  is  established, 
d d m n 


A last  example  concerns  the  Lucas  numbers.  Let  us  show  that 


(5.13) 


,n+l 


L L , - L = (-1)  . 5 

n+1  n-1  n ' ' 


From  (5.1)  and  (5.6)  we  have 

L . . F 


( ."  ) = ( 


, F 
n-1  n 


)(.)• 


n+1 


n+1 


Again  for  = L^  but  from  (5.  4)  with  k = -1  we  get  that 


^n-l\  ! ^n-1  ^n 


( = ( 


n+1 


)(■  ) 


because  L ^ = -1,  L^  = 2.  The  last  two  relations  combined  give 


-1  2 


n n+1  n n+1 


Passing  to  discriminants  and  using  (5.7)  we  obtain  (5.13). 


6.  Some  composite  numbers  that  are  not  Fibonacci  pseudoprimes.  We  have  defined  a 


number  n as  a Fibonacci  pseudoprime  (F.  Psps)  if  it  is  composite  and  satisfies  L^=  1 (mod  n). 


F.  Psps  are  rare:  We  have  seen  that  there  are  only  seven  F.  Psps  < 9161.  It  would  seem  of 

Interest  to  exhibit  some  composite  n which  are  not  F,  Psps.  A modest  beginning  in  this  direc- 
tion are  the  following  results. 

Theorem  4.  The  numbers 

(6.1) 


2 ,(k>l) 


are  not  Fibonacci  pseudoprimes.  Actually 


L^k  = 2*^-1  (mod  2*^) 


-9- 
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(6.3) 


Theorem  5.  If  p Is  an  odd  prime  such  that 


L 1 (mod  p ) , 
P 


L 1^  ^ 1 (mod  p ) ^ k > 1, 


■*  k 

hence  p Is  not  a Fibonacci  pseudoprime. 


For  brevity  we  omit  the  proofs  which  might  be  given  elsewhere.  We  rather  wish  to  dis- 
cuss the  assumption  (6.  3). 

Computer  computations  made  by  George  Logothetis  (November  1976)  show  that 

(6.5)  (mod  n^)  if  2<  n < 7611, 

whether  n is  prime  or  composite.  He  computed  the  remainder  r^,  hence 

2 2 

(6.6)  L = r ( mod  n ) , 0 < r < n , 

n n = n ' 

for  all  n such  that  2 < n < 7611,  with  the  following  results. 

1.  The  remainders  r =0  and  r =1  were  never  found.  This  result  led  us  to  formulate 

n n 

Conjecture  2 of  our  Introduction. 

2.  The  value  = 2 appeared  only  if  n = 0 (mod  24). 

3.  For  n = 24k  he  found  that  = Z precisely  for  the  following  100  values  of  k: 

k=l  2345689  10  12 
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This  is  remarkable  numerical  evidence.  From  generally  large  values,  the  remainder  r^  In  (6.6)  j 

drops  down  to  = 2 for  n = 24k  and  values  of  k as  listed.  We  also  mention  that  the  last  j 

1 

Lucas  number  has  1591  digits.  j 

2 

From  the  identity  - 2 = SIF^^)  [2,  Identity  1^^  on  p.  59]  it  follows  that  = 

5{F^2k^^-  Therefore  - Z = 0 (mod  (24k)^)  if  and  only  if 

(6.7)  Fj2|j  = 0 (mod  24k)  . 

From  the  computer  results  above  we  see  that  (6.7)  holds  for  the  100  values  of  k listed  above, 

and  does  not  hold  for  the  other  values  of  k < [7611/24]  = 317  . | 

i 

i 
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